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Exotic topologically protected zero modes with parafermionic statistics (also called fractionalized 
Majorana modes) have been proposed to emerge in devices fabricated from a fractional quantum 
Hall system and a superconductor. The fractionalized statistics of these modes takes them an 
important step beyond the simplest non-Abelian anyons, Majorana fermions. Building on recent 
advances towards the realization of fractional quantum Hall states of bosonic ultracold atoms, we 
propose a realization of parafermions in a system consisting of Bose-Einstein-condensate trenches 
within a bosonic fractional quantum Hall state. We show that parafermionic zero modes emerge at 
the endpoints of the trenches and give rise to a topologically protected degeneracy. We also discuss 
methods for preparing and detecting these modes. 
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In recent years the concept of topological order has 
revolutionized the way we understand quantum phases of 
matter. Topological phases in one- and two-dimensional 
systems are particularly interesting as the nontrivial ex¬ 
change statistics of particles allows for exotic states of 
matter. For example, Majorana zero modes can emerge 
at boundaries of one-dimensional topological supercon¬ 
ductors [T] or in two-dimensional semiconductor het¬ 
erostructures HSl; see also suggestive experimental sig¬ 
natures in Refs. [SHin]- These topological modes have 
been a subject of intense interest due to their potential 
applications in quantum computation, although they do 
not support universal quantum computing with braiding 
alone [IIHIl]. In two dimensions, certain fractional quan¬ 
tum Hall (FQH) states have been proposed to manifest 
emergent non-Abelian excitations with universal braid¬ 
ing statistics [HHIll. However, an experimental confir¬ 
mation of such emergent non-Abelian anyons has so far 
remained elusive 

Recently, it has been proposed that one can en¬ 
gineer non-Abelian excitations by adding defects in 
the form of ferromagnet-superconductor interfaces at 
the edges of adjacent Abelian FQH states [23H26] . 
The domain wall at their interface binds exotic zero 
modes with parafermionic commutation relations. These 
parafermionic zero modes are associated with a ground- 
state degeneracy that is exponential in the number of do¬ 
main walls (defects). As with any non-Abelian anyon, the 
exchange of the defects binding these parafermionic op¬ 
erators generates a unitary rotation in the ground-state 
subspace. However, the set of operations for quantum 
computation available through such exchanges is richer 
than that available through Majorana exchange [ 23 ]. 
Parafermions may also be realized in bilayer quantum 
hall systems, where the role of the superconducting- (or 
BEC-) induced coupling is played by an interlayer tunnel¬ 



Figure 1: A quasi-one-dimensional finite trench, i.e. a poten¬ 
tial dip, is created by spatially modulating the intensity of 
the laser used to create the dipole trapping potential confin¬ 
ing the atoms to the plane (a cut midway along one trench 
is shown). Within the trench, the two-dimensional density 
of trapped bosons deviates from the FQH filling fraction and 
gives rise to a BEC. 


ing term |28l .'' ]. For a recent proposal on parafermions, 
see also Ref. [29] . 

While existing proposals are based on an experimen¬ 
tally challenging combination of FQH and superconduct¬ 
ing systems of electrons, rapid advances towards creat¬ 
ing a bosonic FQH state open new opportunities to real¬ 
ize topologically nontrivial states in the context of ultra¬ 
cold bosonic systems [301133] . Unlike fermionic systems, 
where a condensed state requires pairing of fermions, 
systems of ultra-cold bosons form Bose-Einstein conden¬ 
sates (BECs) without any additional pairing interaction. 
Thus, if a FQH state is realized in bosons, adding a Bose 
condensed state to such a system can be expected to be 
simpler than in the corresponding fermionic implementa¬ 
tion. 

In this Letter, we propose to realize a bosonic analog 
of a FQH-superconductor system by inserting a trench 


2 


containing a BEC in the middle of a bosonic FQH system, 
as shown in Fig. Such a trench of BEC can be created 
by introducing a potential well that would trap a high 
density of bosons as compared to the FQH region outside. 
We will show that the introduction of such a BEC trench 
induces a novel state in the quantum Hall edge with a 
pairing gap similar to the states with superconductors 
in contact with the FQH system. Further we will show 
that systems containing a pair of such trenches feature a 
ground state degeneracy arising from the topology of the 
underlying quantum Hall system. We explicitly construct 
fractionalized Majorana fermions as localized zero modes 
at the endpoints of the trenches. 

Model .—Let us consider a bosonic FQH state with a 
filling fraction v = 1/m, where m is an even integer. 
While the bulk of the FQH fluid is incompressible, its 
boundaries support gapless edge excitations with frac¬ 
tionalized charge, or boson number in this case, and 
statistics [SH |35] . The edge modes of a FQH state can 
be described by a bosonized field that carries bo¬ 

son number 1/m and satisfies a nontrivial commutation 
relation [5^ 

giv(€)giV3(€') _ g*:^ 

where ^ and are the coordinates along the boundary. 
Equivalently, the chiral field (f satisfies the chiral com¬ 
mutation relations. The creation operator for bosons is 
given by note that bosonic operators commute 

according to the above algebra. Furthermore, the density 
of bosons on the edge is given by p(^) = d^ipj2'K. The 
Hamiltonian iJo = ^ /p {d^(p)‘^, along with the chiral 
commutation ~ ^0 describes a 

free chiral edge mode propagating at velocity v along the 
boundary T of the FQH state. 

As shown in Fig.[^ we consider a BEC trench inside a 
bosonic FQH state that can be introduced by varying the 
potential. The topological phase realized by this FQH 
system will be characterize by a topological degeneracy. 
To create a Hilbert space that allows us to access this 
degeneracy, we will need to consider a system containing 
two such trenches in the FQH system [two copies of the 
trench in Fig. [^. The total Hamiltonian describing the 
boundaries of the FQH state at the edges of the two BEC 
trenches is given by 

H = Hi{ipi^Qi) + H2{f2,Q2) ( 2 ) 

with Qi + Q 2 = 0 mod 1, 

where the subscripts denote the trenches, and Qi is the 
fractional number (“charge”) of quasiparticles modulo 1 
on the edge of the ith trench. While each edge can have a 
fractional charge, the total number of quasiparticles must 
add up to an integer bosonic charge, i.e. Q 1 + Q 2 = 0 mod 
1. We show that, under certain conditions, the Hilbert 
space of this system manifests degenerate ground states 
characterized by parafermionic zero modes at the ends of 


the BEC trenches. The defects at each end of the trench 
thus act as non-Abelian anyons. This is the central result 
of this Letter. We also analyze the robustness of this 
degeneracy to realistic experimental imperfections. 

Single-trench Hamiltonian .—For simplicity, we first fo¬ 
cus on a single edge, with the Hamiltonian Qi), 

without the constraint in Eq. Q, but will impose it later 
in our discussion of the degeneracy of the two-trench sys¬ 
tem. For now, we also drop the subscript 1 for notational 
simplicity. The tunneling between the BEC field and 
the edge states on the trench is described by 

iLtun = -f h.c. (3) 

Note that only a boson, i.e. in quasiparticles bound to¬ 
gether, can directly tunnel to the BEC. The effect of the 
BEC on the edge can be understood by first expanding 
the BEC field as ^'(a;) = 'I'o -f d4'(a;), where i5'I'(a;) is 
the boson fluctuation field in a three-dimensional BEC, 
and is described by a quadratic Hamiltonian. The fluc¬ 
tuations of (5'I'(a;) can be integrated out to obtain an 
effective self-energy for the edge induced by the BEC, 
E = -h E(2), where = -^'qA /pd^ -F h.c. 

is linear in A, while E*^^) is second order in A and is pro¬ 
portional to the Green’s function of the BEC as we shall 
discuss below. 

The lowest order self-energy term E*^^^ does not have 
the form of a cosine term in the sine-Gordon Hamiltonian 
m known to qualitatively modify the chiral edge state 
by opening a gap. The effect of this unconventional term 
can be eliminated by shifting the density p(^) of the edge 
via p(^) —>■ p(^) -l-po for some nonzero po- The density of 
particles at the edge can be tuned by changing the trap 
potential, which in turn changes the chemical potential in 
the BEC. Such a shift can be accommodated by shifting 
the chiral boson field as (p(^) —>■ (p{^) 2 ttpo^. This shift 
introduces a spatial dependence into the lowest order self¬ 
energy E(1) ^ -^-oA ^ g^^ ^nd 

leads to a momentum mismatch |38j . which suppresses 
the effect of the first-order self-energy E^^^. 

Therefore, at a nonzero po, the effect of 
the self-energy is dominated by the second- 
order term E*^^\ which can be written as 
y;(2) _ _^J'^^^(j^'g±jm(<p(?)-l-27rpo€)T*m(v(€')-l-2-n-po5') x 

where a,/3 = 0,1 corresponds to 
the upper/lower sign of the corresponding and 
(i5^'o The correlator in E*-^^ is 

related to the Green’s function Gap of 

the BEC, where r(^) is the position vector of the 
point A more careful discussion of the perturbation 
theory is given in the Supplemental Material [39]. The 
Green’s functions are highly peaked at r(^) « and 
can hence be split into a dominant local part and a 
perturbatively small non-local part. In the remainder 
of the paper, we will first discuss the physics of the 
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Figure 2: The top-view of a BEC trench within a FQH state. 
The trapping potential is engineered such that the boson den¬ 
sity is uniform across the shaded area but vanishes in a small 
region of size I near the endpoints. Quantum fluctuations of 
the BEC couple opposite edges in the shaded region. 


dominant local part where |f(‘C) ~ ^ ^ with W 

the width of the trench, and comment on the effect of 
the perturbative non-local part near the end. 

As shown in Fig. pairs of points directly across 
the trench interact through since they are sepa¬ 
rated by a distance comparable to W. It is convenient 
to define the coordinate x S [0, A] along the trench of 
length L, together with the left- and right-moving fields 
^l{x) = '^{x) and ^Pr{x) = (p{2L — x) propagating along 
the top and bottom edge, respectively. With this def¬ 
inition, the oscillating phase in has a straightfor¬ 
ward interpretation; For example, taking ^ = x and 
^' = 2L — x', the phase corresponding to the Green’s 
function Gqi becomes \ The presence of the lat¬ 

ter term ensures momentum conservation [38j ; only such 
combinations contribute significantly to Physically 

the only terms that satisfy these constraints are nor¬ 
mal particle-particle scattering terms on the same edge 
or anomalous particle-hole scattering between opposite 
edges. The former simply renormalizes the velocity, while 
the latter is analogous to superconducting pairing. Tak¬ 
ing the above constraints into account, the dominant, 
and local, contribution to the effective interaction is 

/ L-l 

dx cos{2m(j)), (4) 

where we have defined the (non-chiral) fields (j^ix) 
and 9{x) by <fL/R{x) = (f’ix) T S{x), and the cou¬ 
pling coefficient A cx A^. The new non-chiral 
fields are self-commuting at all points, but satisfy 
[(j){x),6{x')] = i{'K/m)Q{x — x'). 

The term cos{2m(f>) represents the effective pairing be¬ 
tween the edges and is reminiscent of the sine-Gordon 
model m- The resulting Hamiltonian that describes the 
BEG trench (Fig.[^ can be written as 

HeS ~ ^ y dx[{d^(l))‘^ + ^ J dxcos{2m(j)), 

( 5 ) 

where the continuity of dxf{x) near x ^ 0,L transforms 
into the boundary condition dx4>{x = 0, L) = 0. We 
note that, in addition, we are required to preserve the 


boundary conditions for 0 at the endpoints: 0(x = L) = 0 
follows from the definition, while 9(x = 0) = ttA/to, 
where N is the total number of quasiparticles. The latter 
is due to the fact that the total density of bosons on 
both edges is p = —dxdjTT. However, similar to other 
restrictions on the Hilbert space, we will restore these 
boundary conditions at the end of the calculation. At 
large A > 0, the sine-Gordon model [Eq. ([^] supports 
several ground states \p) characterized by the expectation 
values 

p = 0,l,--- ,2m-l, (6) 

for X away from the edges and in the limit of large L[37]. 

To restore the appropriate boundary conditions for 
9{x = 0, L), we notice that 9{x = L) commutes with (j){x) 
and Heff, hence it can be set to any value [9(x = L) = 0 
in our case] without consequence for the single trench. 
On the other hand, 9(x = 0) = irN/m obeys a non-trivial 
commutation relation with ^i(x) such that [N, (j)\ = —i\ 
however, R = commutes with Heff- The oper¬ 

ator R, despite being a symmetry of the Hamiltonian, 
transforms (/>(x) in a non-trivial way as 

R(j)(x)R^ = Mx) H-. (7) 

m 

It then follows that, including the boundary conditions, 
different values of (j) in Eq. § indeed correspond to the 
same energy. Exploiting the above symmetry, we can 
describe the ground states in a basis that also makes the 
operator R diagonal as 

2m—1 

\n}= Y, \P ). ( 8 ) 

p—0 


which satisfy R\n) = . Physically, n corre¬ 

sponds to the number of quasiparticles modulo pairs of 
bosons on the edge of a given trench. In our model thus 
far, this is a well-defined quantum number since the co¬ 
sine term in Eq. ([^ transfers only pairs of bosons from 
the condensate. Note that N takes 2m distinct eigenval¬ 
ues. This allows us to associate a degeneracy of y/2m 
with each of the endpoints of the BEC trench (the ‘de¬ 
fects’) in the limit of a large number of trenches where 
restrictions on the Hilbert space may be safely ignored. 
We will discuss Hilbert space constraints for a small num¬ 
ber of trenches in the next section. Finally we remark 
that the fractional part of the boson number on a trench, 
which is invariant under the addition of single bosons, is 
Q = N/m mod 1, while Qt = [N/m — Q] S {0,1} defines 
the boson parity. The '/2m degeneracy thus also requires 
protection by a Z 2 symmetry due to boson parity. 

Degeneracy .—So far we have focused on the spectrum 
of a single trench without the physical constraints of 
the Hilbert space. With a finite nnmber of trenches, 
the total fractional charge Q should be 0 mod 1 |45j . 
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Considering the double trench model of Eq. ([^, the 
fractional boson numbers on the two trenches satisfy 
Qi = —Q 2 = 0, 1/m, ■ ■ ■ ,1— 1/m, while their boson par¬ 
ities Qb i = 0,1 and Qt ,2 = 0 ) 1 are unconstrained, which 
yields a total degeneracy ol D = Am for two trenches. 
More generally, a system of k trenches with boson-parity 
conservation intact has the degeneracy 


D^2{2m)^-\ (9) 

We shall discuss the (topological) robustness of this 
degeneracy, which partially survives the Z 2 symmetry 
breaking, after we explicitly construct operators that 
span the 4m degeneracy present in the two trench case. 

Parafermion operators. —We now construct the oper¬ 
ators spanning the above two-trench degeneracy. As re¬ 
marked earlier, this degeneracy is spanned by the ex¬ 
change of quasiparticles between the endpoints of the 
trenches. Let Ui {U/) be the exchange operator between 
the left (right) endpoints of the two trenches. They can 
be expressed as 


Ulr ^ ( 10 ) 

where and represent quasiparticle vertex operators 
acting on the top trench and bottom trench, respectively, 
and projected onto the ground-state sector. For exam¬ 
ple, t/j, = adds a quasiparticle on one of the 

two ends of the first trench with the projection onto the 
ground state implicit. The operator B on the second 
trench can be defined similarly by ipi —> ip 2 ', however, 
we must insure that the exchange of quasiparticles be¬ 
tween the two trenches respects the exchange statistics. 
This can be done consistently by including a Klein factor 
as b\^ ^ where A^i is the number of 

quasiparticles on the first trench. With this construction, 
we now focus on the operator T defined above. Note that 
the chiral field at the ends is given by tpi = (f>i — di. As 
we discussed earlier, 9i{L) — 0 and 0i(O) = nNg^i/m 
mod 27r with Ng^i = m(Qt,p -I- Qi) the total number of 
quasiparticles modulo 2m on the first trench. Since we 
are interested in the ground state sector, the field is 
roughly assumed to be pinned according to Eq. Q over 
the entire edge as (/'i(O) « '/’i(L) « 27rfV0p/m with inte¬ 
ger N^^i. We find 


±l — e , — e , 


which satisfy the algebra 


{tIY^ 


= 1 , 


r/Tt 


= e 


T/r/ 


( 11 ) 


( 12 ) 


These relations describe parafermions, a generalization 
of the fermionic algebra mmiiD], see also Ref. m- 
The operators Bi ^ can be defined similar to Eq. 0 
by including the Klein factor exp{iTrNg^i/m) explained 
above, and with Ngj —)■ Ng^2 = m{—Qi + Qb,i) and 


-^ 0 . 2 - The parafermionic algebra implies that 
(^t) 2 m, _ which yields 2 m degenerate ground states in 
a sector with a fixed total parity [the operator U does not 
change the total parity according to Eq. @]. With the 
two-fold degeneracy due to the total parity, one recovers 
the full 4m degeneracy of the system. 

Robustness. —Heretofore, we have focused on the 
Hamiltonian in Eq. ([^. In principle, however, a sin¬ 
gle boson can tunnel to or from the BEC, which might 
arise from terms of the type V = cos[m(j){x = 0, L)] near 
the endpoints of the trench. In fact, introducing the per¬ 
turbation V within first order degenerate perturbation 
theory reduces the 4m fold degeneracy to D' = m/2-, see 
the Supplemental Material [39] for details. While such 
a term clearly breaks the boson parity symmetry (hence 
(Qf).i: Qb. 2 ) are no longer good quantum numbers), it also 
breaks a two-fold degeneracy of the fractional number of 
quasiparticles Qi. The remaining m/2-fold degeneracy is 
topologically protected, and we find in general that the 
degeneracy for k trenches is given by 


D' = {m/2) 


fc-i 


(13) 


This gives a quantum dimension of sJm/2 for the de¬ 
fects ending the BEC trenches. However, if the vacuum 
regions (Fig. are sufficiently large, single-boson tun¬ 
neling is suppressed, and one recovers the degeneracy in 

Eq. Q. 

We also briefly remark that the long-range fluctuations 
of the BEC can be considered effectively as long-range 
tunneling, and can provide another mechanism to break 
the boson-parity symmetry, while the quantum dimen¬ 


sion in Eq. (13) will not be affected. 


Preparation and detection. —There is some evidence 
that & V = 1/2 fractional Chern insulator has a continu¬ 
ous transition to a BEC, which allows for quasi-adiabatic 
preparation of the former [42] . A similar procedure may 
exist for FQH states. The parafermions in the FQH state 
can be prepared by starting with a small island which can 
be grown to a trench in linear time by modulating the 
laser beams. Furthermore, Bragg spectroscopy can pro¬ 
vide direct information about the topological phase of 
the system. For example, one should observe a zero bias 
peak at the endpoints of the trench [33]. They can also 
be probed using braiding, which corresponds to the topo¬ 
logically protected manipulation of the underlying quan¬ 
tum information and which requires dynamically chang¬ 
ing the geometry of the system, bringing different sets 
of parafermionic edge modes in close proximity to each 
other [331 mi- Such dynamical changes can easily be 
achieved by dynamically changing the laser beam used 
to create the BEC trenches. 

It is worth pointing out that a theoretically simpler but 
experimentally more challenging approach would proceed 
by analogy with Ref. [33]: A BEC of diatomic molecules 
can be coherently dissociated into pairs of atoms, which 
readily gives the effective Hamiltonian in Eq. Q. 
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Conclusion .—In this work, we have considered a BEC 
trench in a FQH liquid, and showed that, in a certain 
regime, the combined system is in a topological phase, 
which is identified by the zero mode operators at the 
endpoints of the trench. These zero modes are shown to 
be parafermions, a generalization of the usual fermionic 
or bosonic algebra. We have also derived the topological 
degeneracy of the parafermionic modes, and examined 
their robustness against local perturbations. 
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Supplemental Material 


SI. EFFECTIVE INTERACTION TO THE 
SECOND ORDER 

We start with the interaction in Eq. (3). It is more con¬ 
venient to work in Euclidean space where the interaction 
takes the form 

- A y 0 + C.C. (SI) 

We expand the condensate held as = dtQ-l-Jdt similar to 
the main text. The effective interaction, after integrating 
out the BEC held, is formally given by 

g-e[¥)(r,{)] _g'I'oA JdTd{e™‘^-|-c.c. 

X /gA \ 

\ /(S'!/ 

where the average {■)s'}> is taken over the huctuations of 
(Jtlt. To hrst order, the effective interaction is given by 
the prefactor in the above equation. To second order, the 
effective interaction is 

r 

^ / drclr'd^d^' x 

^^m{v{r,i)+v{r',i'))^^2^mpo{i+e)Q^^ T, t', r'] -f C.C. 
g^n^(<,.(r.^)-v(r^^'))e*2™po(^-^')G<gp[^^ r', r'] -f c.c.j, 

(S3) 

where r = r(^) and r' = r{^'), and 

Goo(T,r,0,0) = (<5vI/(T,e)5^^(0,0)), 
Goi(r,r,0,0) = {S'S{t,OS^{0,0)), (S4) 

are the normal and anomalous Green’s functions of the 
BEC, respectively. While these functions are nonlocal 
both in space and (imaginary) time, they are highly 
peaked around (t, r) —> 0, and thus the interaction is as¬ 
sumed to be dominantly local plus small non-local terms. 
In the second-order perturbation theory treatment of the 
main text, we have neglected the time-dependence of the 
Green’s functions to simplify the presentation. With the 
interaction approximated to be local in time as well as 
space, one can write the effective interaction as 

0[(p(t,a;)] = J dTT,[ip{T,x)] = j drdxVeS {(p{t,x)), 

. 

where x is the position along the trench defined in the 
main text, and Veg is the effective Hamiltonian, which, 
after dropping the momentum-mismatched terms, be¬ 
comes Eq. (4) in the main text. 


S2. DERIVATION OF THE QUANTUM 
DIMENSION IN THE ABSENCE OF 
BOSON-PARITY SYMMETRY 

In this Section, we perform first-order perturbation 
theory in the absence of boson-parity symmetry and de¬ 
rive the quantum dimension given in Eq. (13). 

We begin by consider a single trench without any con¬ 
straints, and only later implement the constraint on the 
physical Hilbert space. The parafermion operators inro- 
duced in the main text allow us to define a state with n 
quasiparticles on the edge of the trench as 

|n) = (Tt)”|0), (S6) 

where |0) is the state with no fractional charge. The 
single-boson operators on the two ends of the trench are 
defined as 

bi/r={Tl^y\0), (S7) 

at the corresponding left/right endpoints. We have 

br\n) = \n + m), (S8) 

by definition, and 

k\n) = bi (rt)”|0) = (-!)" (Tt)”6;|0) 

= cr(—l)"|n-I-m), (S9) 

where a = {m\bi\0) = (0|&);&;|0); in deriving the 

above equation, we used the commutation relation of 
parafermion operators in Eq. (12) of the main text, and 
specifically biT} = —T}bi. We also note that the quan¬ 
tum numbers in brackets should be understood modulo 
2m. 

Using the definitions given in the main text, one can 
see that the operators bi and b^ are Hermitian, and 
a = (—1)™/^ is a real number (m is even). Therefore, 
a perturabtion of the Hamiltonian that adds a single bo¬ 
son, and thus breaks the boson parity, has the form 

H' = abi + j3br, (SIO) 

with a and /? being arbitrary real numbers. This Hamil¬ 
tonian mixes the states \n) and |n-|-m) as {n+m\H'\n) = 
(—1)”ctq; -I- (3. Thus, the eigen-energies are split as 

±[(-l)Va + /3], 

which implies that the number of minimum-energy states 
is reduced by a factor of 4, from 2m (see the main text) 
to m/2. 

In the next step, we consider two trenches, and impose 
the physical constraint that the total fractional charge is 
trivial, i.e., = 0 mod m where the superscripts 

denote the two trenches. In particular, this implies that 
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(_!)«( ' _ (— 1 )”* \ A similar perturbative treatment 
yields the energy splitting 


± 





(Sll) 


with and the corresponding perturbation param¬ 
eters for each trench. The splitting of energy reduces the 
total degeneracy down to m/2. (Note that there is no 


freedom in choosing the relative parity of and 
due to the mixing of and -I- m in the true ground 
states of the perturbed Hamiltonian.) A straightforward 
generalization to arbitrary number of trenches shows that 
the quantum dimension of the endpoint defects in the ab¬ 
sence of boson-parity symmetry is given by Eq. (13) of 
the main text. 



